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We report the result of a search for the pair production of the lightest supersymmetric partner of 
the top quark (ti) in pp collisions at a center-of-mass energy of 1.96 TeV at the Fermilab Tevatron 
collider corresponding to an integrated luminosity of 5.4 fb _1 . The scalar top quarks are assumed to 
decay into a b quark, a charged lepton, and a scalar neutrino (u), and the search is performed in the 
electron plus muon final state. No significant excess of events above the standard model prediction 
is detected, and improved exclusion limits at the 95% C.L. are set in the the (Mj l ,Mp) mass plane. 
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Supersymmetric theories [l| predict the existence of 
scalar partners for each of the standard model (SM) 
fermions. In the minimal supersymmetric standard 
model (MSSM) 0, the mixing between the chiral states 
of the scalar partners of the SM fermions is greatest for 
the partners of the top quark due to its large Yukawa cou- 
pling [3\. Thus, it is possible that the scalar top quark 
(ti ) is the lightest squark and has the largest production 
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cross section. If i?-parity [l| is conserved, then scalar top 
quarks would be produced by pp collisions in pairs with 
the dominant processes being quark- antiquark annihila- 
tion and gluon fusion Q. 

In_this letter we report on a search for the production 
of titi pairs in the bbe^y^vv final state. We assume that 
the t\ has a 100% branching fraction in this three-body 
decay mode with equal fraction to each lepton type, that 
i?-parity is conserved, and that the sneutrino (y) is the 
lightest supersymmetric particle or decays invisibly into 
a neutrino and a neutralino (x?)- This analysis uses data 
corresponding to an integrated luminosity of 5.4 fb _1 col- 
lected using the DO detector operating at the Fermilab 
Tevatron collider at y/s = 1.96 TeV. The data were col- 
lected from April 2002 through June 2009. The DO Col- 
laboration has previously searched for top squark 
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pair production in the final states bbl^l^vv where the 
lepton pair is ee, /j,/x, or e[i. Two of these earlier searches 
used subsets of this data set corresponding to integrated 
luminosities of 0.43 fb _1 and 1.1 fb -1 , while the earliest 
search used data from the Tevatron Run I, corresponding 
to an integrated luminosity of 0.11 fb _1 . Searches for top 
squark pair production in the bbi^i^vv final states have 
also been reported by the CDF collaboration Q and by 
the ALEPH, L3, and OPAL Collaborations B. 

The main components of the DO detector [9| include 
a central tracking system located inside a 2 T super- 
conducting solenoid. The inner-most tracking element 
is the silicon microstrip tracker (SMT), followed by a 
scintillating fiber tracker. These two detectors together 
measure the momenta of charged particles. The tracking 
system provides full coverage in the azimuthal (0) direc- 
tion for 1 77 1 < 2, where the pseudorapidity rj is defined as 
7] = — ln(tan#/2) and 9 is the polar angle with respect to 
the proton beam direction. Outside the solenoid is the 
uranium/liquid argon calorimeter which is divided into a 
central calorimeter and two end-cap calorimeters. Each 
of these three calorimeters have electromagnetic layers 
followed by hadronic layers. The outermost component 
of the detector is the muon system, which consists of 
proportional drift tubes and scintillator trigger counters, 
followed by 1.8 T iron toroids and two additional layers 
of drift tubes and scintillators. Events are selected for of- 
fline analysis by a three-level trigger system. All events 
are required to pass one of a suite of single-electron trig- 
gers or single-muon triggers using information from the 
tracking system, the calorimeter, and the muon system. 

For each event, a primary pp interaction vertex is de- 
fined. If more than one vertex is reconstructed, the pri- 
mary vertex is taken to be the vertex least consistent 
with originating from a soft collision. The location of 
the primary vertex along the beam direction is required 
to be within ±60 cm of the detector center. 

Jets are reconstructed using the DO Run II 
cone algorithm [l(| with cone of radius 1Z = 
J (A0) 2 + (Ay) 2 < 0.5, where y is the rapidity. 
Jet energies are calibrated using the standard DO pro- 
cedure [llj. A jet is retained in an event if it has 
transverse energy Et > 20 GeV, \rj\ < 2.5, and if 
Aft(jet,electron) = ^(A^) 2 + (Arj) 2 > 0.5. No require- 
ment on the number of jets is applied. 

Electrons are required to have transverse momen- 
tum pj, > 15 GeV and \rj\ < 1.1. They are re- 
quired to be isolated, defined as having [E tot (QA) — 
E em {0.2)]/E E m{0-2) < 0.15 where E tot (QA) is the to- 
tal calorimeter energy in a cone of radius TZ — 0.4 and 
Eem(0-2) is the electromagnetic energy in a cone of ra- 
dius TZ = 0.2. In addition, the shower development in 
the calorimeter is required to be consistent with that of 
an electromagnetic shower both transversely and longi- 
tudinally. An eight-variable likelihood function is con- 
structed to further distinguish between electromagnetic 



and hadronic showers. The output of this function ranges 
from to 1, and electrons are required to have a likeli- 
hood value greater than 0.85. Electromagnetic showers 
associated with electrons are also required to match a 
central track within A77 < 0.05 and Acf> < 0.05 of the 
electromagnetic cluster. 

Muons are required to have transverse momentum 
Pj, > 10 GeV and \r)\ < 2. They are required to have both 
drift tube and scintillator hits in the muon system and to 
match a track in the central tracker. If the central track 
includes hits in the SMT, the distance of closest approach 
(DCA) between the muon track and the primary vertex 
is required to be less than 0.02 cm. If there are no SMT 
hits, then the DCA is required to be less than 0.2 cm. 
Muons are also required to satisfy isolation requirements 
in both the calorimeter and the central tracker. For the 
calorimeter isolation, the transverse energy in the cone 
TZ < 0.5 around the muon track divided by p^ must be 
less than 0.15. For the central tracker isolation, the sum 
of the transverse energy of the tracks in the hollow cone 
0.1 < 1Z < 0.5 divided by pj, must be less than 0.15. 

Events are required to have exactly one electron and 
one muon with opposite charge and to have a min- 
imum separation between the electron and the muon 
AlZ(e, n) > 0.5. The missing transverse energy {$ T ) 
is calculated from the calorimeter energy corrected for 
the jet and electron calibrations. It is then adjusted to 
account for the transverse momentum of the muon. All 
retained events are required to have $ T > 7 GeV. We re- 
fer to this preliminary set of selection requirements as 
the preselection. 

Signal Monte Carlo (MC) events are generated in a 2-D 
grid, i.e., for ii masses ranging from 100 GeV to 240 GeV, 
and for v masses ranging from 40 GeV to 140 GeV, each 
in 10 GeV steps. For each point, the MSSM decay param- 
eters are calculated with SUSPECT [H| and SDECAY [HI]. 
madgraph/madevent [lj] is used to gene rate four- 
vectors for the signal events with pythia [15] providing 
the showering and hadronization. The next-to-leading 
order (NLO) cross section for t\ pair production is cal- 
culated by prospino 2.0 [B] with the CTEQ6.1M [T3] 
parton distribution functions (PDFs). The calculations 
are performed with the factorization and renormalization 
scales set to one, one half, and two times the t\ mass to 
determine the nominal value and the negative and posi- 
tive uncertainties. The scale factor uncertainties are com- 



bined quadratically with the PDF uncertainties [12, |18[ 
to give the total theoretical uncertainties for the signal 
cross sections. 

The dominant SM backgrounds for this decay are 
Z/7* — > tt with t — > lv; diboson production includ- 
ing WW, WZ, and ZZ; top quark pairs; W + jets; and 
instrumental background coming from multijet (MJ) pro- 
cesses where jets are misidentified as electrons or contain 
muons that pass the isolation criterion and with $ T aris- 
ing from energy mismeasurement. All the background 
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processes in this analysis except for MJ are modeled us- 
ing MC simulation. Vector boson pair production is sim- 
ulated with pythia, while all other backgrounds are sim- 
ulated at the parton level with ALPGEN [19| , with PYTHIA 
used for hadronization and showering. In order to sim- 
ulate detector noise and multiple pp interaction effects, 
each MC event is overlayed with a data event from ran- 
domly chosen pp crossings. 

MC correction factors determined from data are ap- 
plied to make distributions consistent between data and 
MC. These corrections include factors for the luminosity 
profile, beam spot position, muon and electron identifi- 
cation efficiencies, boson transverse momentum, and jet, 
electron, and muon energy resolutions. 

The MJ background is estimated from a selection of 
data events not overlapping with the search sample and 
is selected by inverting the electron likelihood and muon 
isolation requirements. This sample is used to determine 
the shape of the MJ background. Because most same- 
sign di-lepton events come from MJ processes, we obtain 
the normalization factor by taking the ratio of the num- 
ber of same-sign events that pass the likelihood and iso- 
lation requirements to the number of same-sign events 
that fail these requirements. To remove W+jet events 
from the MJ same-sign sample, we make the additional 
requirement $ T < 20 GeV, since W^+jets events tend to 
have large $ T . We also correct this ratio for non-MJ SM 
processes that produce like-sign leptons, using the MC 
samples. 

Data events are required to satisfy at least one of a 
suite of single-electron or single-muon triggers. The effi- 
ciency of the combination of the single-electron triggers is 
measured using a subset of the search sample in which at 
least one of the single-muon triggers fired, and vice- versa 
for the single muon triggers. The combination of these 
two efficiencies, taken to be the overall trigger efficiency, 
is then applied as a correction to the MC samples. 
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FIG. 1: (color online) $ T versus A<f>(e,fi) for Z/-y* -» tt MC 
events. Events inside the black box in the lower right corner 
are removed. The sizes of the boxes are proportional to the 
number of events in the underlying cell. 



The mass difference, AM = Mj — Mp determines 
the kinematics of the final state. A larger AM will 
lead, on average, to larger $ T , larger jet energy, and 
higher px charged leptons. We divide the range of 
AM into a "large-AM" region (AM > 60 GeV) 
and "small-AM" region (AM <60 GeV). To illustrate 
these regions, we have chosen two benchmark points, 
(M fl , M>) = (200 GeV, 100 GeV) and (110 GeV, 
90 GeV), which will be referred to as the large-AM 
and small-AM benchmarks, respectively. Since there are 
many signal points and their characteristics differ sig- 
nificantly, the analysis strategy is to optimize the signal 
selection as a function of AM. 

For all values of AM, the largest background after pre- 
selection is Z/j* —> tt. A two-dimensional plot of the az- 
imuthal angle between the electron and muon, A<^(e, /i), 
vs. f T for Z/j* — > tt MC events is shown in Fig. [TJ The 
two leptons from Z/j* ~ > tt tend to be back-to-back in 
(j>, and tend to have low $ T . We therefore reject events 
in which A(j)(e, fj) > 2.8 and $ T < 20 GeV and label this 
as "Selection 1". 

Figure compares $J T , electron px , and muon px of 
the data and the sum of all backgrounds at this stage of 
the analysis. The agreement confirms our understanding 
of the SM backgrounds, of the trigger efficiency, and of 
other MC corrections. After selection 1, the three largest 
backgrounds are Z/j* — > tt, WW, and it production. 
To discriminate these backgrounds from signal we create 
for each of them a composite discriminant variable from 
a linear combination of kinematic quantities. We use 
the R software package j2(3] to calculate the maximum 
likelihood coefficients for a generalized linear model 
(GLM) [HI of the form 



5 A = In ■ 



(i 



1 



/' 



A) + £ • x 



(i) 



to discriminate between signal and a specific background 
source A. Here, \i is the probability that an event is 
signal, (3q is a constant, j3 is the vector of coefficients, 
and X is the vector of event kinematic variables. By 
construction, 5 A = when \i = 0.5, and signal-like events 
have positive 5A. The discriminant SZ is constructed 
to separate signal from Z/j* — > tt background, using 
an equal number of signal and Z/j* — > tt MC events 
to determine the coefficients (3o and f3. For X we use 
the following variables: ln($ T ), ln(p^), ln(p T ), A</>(e,/i), 
A(j){e,$ T ), A<j){ii,$ T ), and A(j)(e,$ T ) x A0(/i, $ T ). For 
each value of AM, ranging from 20 to 200 GeV, we use 
the same variables with re-optimized coefficients. We use 
a similar method for creating the discriminants SWW 
and Stt to separate signal from WW and tt backgrounds. 
For SWW we use the variables ln($ T ), ln(pj.), ln(pf,), 
number of jets, A(f>(e,fj,), and ln(WWt a g)- Here WW ta , e 
is the magnitude of the vector sum of p 1 ^, p 1 ^, and $ T , 
which should be close to zero for WW events. For Stt, we 
use the variables ln($ T ), m(p^), ln(pf,), ln(l + Ht), the 
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FIG. 2: (color online) Distribution of (a) $ T , (b) electron 
Pt, and (c) muon pr comparing data and all background pro- 
cesses after Selection 1. The thick dashed and thick solid lines 
represent the small-AA/ and large-AA/ signal benchmarks, 
respectively. 



energy of the second most energetic jet, and WWt ag - The 
variable Ht is the scalar sum of the transverse energies 
of all jets in an event. 

We first apply a requirement using the most effective 
discriminator of the three. For AM < 60 GeV, we re- 
quire Stt > 0. The efficiency of this requirement is 0.95 
for the small-AM signal benchmark and 0.03 for ti. For 
AM > 60 GeV, we require SZ > 0. The efficiency of this 
requirement is 0.96 for the large-AAf signal benchmark 
and 0.01 for Z/"f* —¥ tt. After making these require- 



ments on one variable, we build 2-D distributions of the 
two remaining discriminants. Figure [3] shows these dis- 
tributions for the small-AM benchmark signal and the 
two most significant remaining backgrounds. In calcu- 
lating the signal exclusion confidence limits, we use only 
the bins in the upper right quadrant where the signal is 
concentrated. 
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FIG. 3: (color online) Distribution of SZ versus SWW for (a) 
the small-AM signal benchmark, (Mj^Mp) = (110 GeV, 90 
GeV), MC events, (b) Z/'y* -»■ tt MC events, and (c) WW 
MC events. The top right quadrant is used in the limit-setting 
procedure. 



Tabic U summarizes the expected backgrounds, the ex- 
pected signal, the observed data events, and the selection 
efficiencies. 

The theoretical uncertainties on the signal cross section 
are approximately 20% as discussed above and are the 
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Preselection Selection 1 Selection 2 



Sample 


Events 


Events 


Stt > 
Events 


5Z > 
Events 


Z — T TT 


lOlD ± 1DU 


coo _l_ fil 
DoZ ± 01 


Kl c _1_ KA 


170 1 oil 


Z — T HfJ, 


QQ 1 1 A 1 

66.1 ± 4. i 


ZZ.y ± 6. 1 


ICQ lOfl 


.0 ± 1 . z 


Z — t CC 


OQ O-I-QQ 

zo.z it o.y 


1D.D it O.U 


i n o -t- o o 


U . 1 ± Z. o 


vv z 


19 7-1-1 
1Z. / It 1.0 


10 0-1-1 ^ 
1Z.U it 1.0 


D.o =t U.o 


Q Q _1_ 1 O 

y.o ztz i.z 


T T r T T r 
W W 


z9o ± 62 


208 ± 30 


157 ± 18 


2,61 ± zo 


zz 


2.2 ±0.3 


2.0 ±0.3 


1.0±0.15 


1.1 ±0.16 


ft 


206 ± 28 


204 ± 28 


6.6 ±0.9 


179 ± 24 


w 


70 ±9.2 


67.5 ±9.0 


55 ±7.7 


53 ± 7.4 


MJ 


33 ±9.2 


19.7 ±5.5 


18.4 ±5.1 


1.3 ±0.35 


Background total 
Data 


2191 ±160 
2168 


1195 ±73 
1147 


785 ± 57 
776 


513 ±37 
472 


Small-AM Benchmark 
(110 GcV,90 GcV) 


35 ±5.6 


25.5 ±4.2 


23.8 ±3.9 




Large-AM Benchmark 
(200 GcV,100 GcV) 


55 ±9.3 


53.4 ±9.0 




51.8±8.7 



TABLE I: Expected numbers of background and signal events, and the number of events observed in the data at each stage 
of the analysis. The errors include statistical and systematic uncertainties. For AM < 60 GeV, Selection 2 is Sti(AM — 20 
GeV)> 0, and for AM > 60 GeV, Selection 2 is 8Z{AM = 100 GeV)> 0. 



dominant uncertainties in this analysis. The uncertainty 
on the integrated luminosity is 6.1%. Other systematic 
uncertainties included in the limit setting calculations 
are the lepton identification and track matching efficien- 
cies (5%), the MJ background (27%) scale factor, the 
jet energy calibration (1-2)%, and the production cross 
section uncertainties on all the SM background processes 
(3-10)%. All uncertainties except for those on the MJ 
background and the SM production cross sections are 
treated as fully correlated. All systematic uncertainties 
are included in the limit calculations with Gaussian dis- 
tributions [Hf. 

For AM < 60 GeV, we use the two dimensional his- 
tograms of the positive values of 6Z and SWW in the 
limit setting procedure. For AM > 60 GeV, we use the 
positive values of Sti and SWW. A modified frequentist 
approach [23j is used to determine the 95% C.L. exclu- 
sion limits on scalar top quark production as a function 
of the v and t\ masses, as shown in Fig. HI Also shown 
are the exclusion regions from the CERN LEP experi- 
ments 8]. previous DO searches [5. 6]. and a CDF search 
0. 

In conclusion, we set 95% C.L. exclusion limits on the 
cross section for scalar top quark pair production assum- 
ing a 100% branching fraction to bbl^l^vv using 5.4 fb _1 
of integrated luminosity from the DO experiment at the 
Fcrmilab Tevatron collider. We have excluded stop pair 
production for M~ ti < 210 GcV when M v < 110 GcV 
and the difference M i± - M > 30 GeV. This extends 
the previous limits on the top squark mass by more than 



40 GeV for sneutrino masses less than 90 GeV and the 
limits on the sneutrino mass by more than 30 GeV for 
top squark mass equal to 150 GeV. 
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